Carrier drift velocity and edge magnetoplasmons in graphene 
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We investigate electron dynamics at the graphene edge by studying the propagation of collective 
edge magnetoplasmon (EMP) excitations. By timing the travel of narrow wave-packets on picosec- 
ond time scales around exfoliated samples, we find chiral propagation with low attenuation at a 
velocity which is quantized on Hall plateaus. We extract the carrier drift contribution from the 
EMP propagation and find it to be slightly less than the Fermi velocity, as expected for an abrupt 
edge. We also extract the characteristic length for Coulomb interaction at the edge and find it to 
be smaller than for soft, depletion edge systems. 
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The Quantum Hall Effect (QHE) occurring in 2- 
Dimensional Electron Systems (2DES) relies vitally on 
edges, both to accumulate charge and to support non- 
dissipative chiral currents. The time dependent Hall ef- 
fect takes the form of Edge Magneto-Plasmons (EMP), 
quasi one dimensional gapless elementary excitations 
which are split off from the bulk magneto-plasmon modes 
by the sample boundary (H-[l3j . Ever more closely con- 
fined to the edge as frequency and wavevector are in- 
creased, they are a tool of choice to investigate edge 
structure and dynamics. Seen in both classical]!, E| and 
quantum inertial 2DES [8MT3I . they could also exist for 
electrons in graphene despite the very different gapless 
Dirac dynamics [LJ, [l5|. We present experiments which 
establish firstly that EMP do indeed exist in graphene 
and then, from the propagation properties, we extract 
both the velocity of the chiral current on the edge and 
the characteristic length for edge Coulomb interaction. 

Electrons in graphene obey two-dimensional massless 
rclativistic-type dynamics with speed vp ~ 10 8 cm s _1 
(Fermi velocity). The honeycomb lattice structure adds 
pseudospin-orbit coupling. Energy dispersion is linear 
and gapless around the band crossing "Dirac"_points sit- 
uated at the Fermi level for charge neutrality [16j. Charg- 
ing sweeps the Fermi level through the crossing to make a 
smooth transition from electron to hole behaviour. The 
electronic edge is then never defined by depletion (the 
"soft" edge of the usual gapped 2DES) and the electrons 
should always feel the abrupt work function potential of 
the sample edge. On approaching the boundary in per- 
pendicular magnetic field the energies of the edge Landau 
states increase because of confinement and eventually 
emerge above the Fermi level. The resulting confinement 
force is balanced by carrier drift along the edge which 
constitutes a chiral current with vanishing backscatter- 
ing, the essential ingredient for a quantum Hall effect 



which, in graphene, because of the unusual degeneracies 
induced by the pseudo-spin coupling, displays an anoma- 
lous integer filling factor sequence v = 2, 6, 10, . . . [I?], EH , 
v being the number of electrons per flux quantum. 

Classically EMP appear as propagating periodic vari- 
ations of charge excess/deficit on the edge which, by cre- 
ating an electric field, induce drift currents which rede- 
posit the charge to advance the wave along the boundary. 
Because the charge accumulation is proportional to the 
Hall drift current, the phase velocity is determined by the 
Hall conductivity a xy and electrostatics. On the edge the 
charge drifts at velocity vd, pulling the electrostatic field 
with it and thereby adding Dn to the Hall conductivity 
propagation. In quantum language EMP become a prop- 
agating local variation of chemical potential (filling of the 
Landau level as it emerges above the Fermi level at the 
edge) which drives a local Hall current which by con- 
tinuity is accompanied by charge propagation 19]. The 
inward force driving the electron drift is set by the vari- 
ation of chemical potential with charge imbalance and 
consists of two parts: the slope of the emerging Landau 
level and the electrostatic field from the imbalance which 
give rise respecti vely to the vd and a xy terms in the dis- 
persion relation 19J, l20j 
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where e g- is the effective dielectric substrate screening 
constant, q the wave vector, w a lower cut off length in 
the electrostatics which can be expected to be of the or- 
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der of the magnetic length £3 = (hc/eB) 1 , but could 
also result from lateral spread of the charge imbalance. 
The constant C = 1 for a step function density profile, 
but in general it is a function of the equilibrium pro- 
file at the edgejfjj]. The basic form of Eq. (1) is the 
same for graphene as for the usual 2DES because Lorentz 



force drift motion depends ultimately only on the Lorentz 
transformation to the zero electric force frame for the 
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electrons |2l|. The slope of the Landau level energy, and 
therefore the drift velocity, is, however, dynamics depen- 
dent and scales with the respective Fermi velocities. The 
density profile is determined principally by electrostatics 
with narrow incompressible strips where the density situ- 
ates the Fermi level in the Landau level gaps[22j . Electro- 
statics for a charged equipotential sheet imposes the lim- 
iting form of the edge density profile to be n s (x) ~ x~ 1 ^ 2 
where x denotes distance from the edge, as opposed to 
n s {x) ~ x 1 / 2 for a depletion edge. The repulsive force re- 
sponsible for the drift velocity results from the increase in 
energy W n ~ pi)p ~ *JnKv f/£b of the Landau function 
of level index n = (y — 2)/4 as it is squeezed against the 
hard edge: in graphene, the effective confining electric 

field is Ef = -(l/e)6W n /5x hv F /e£%, where £ B 

characterises the width of the wavefunction. Electrons 
drift along the edge in this potential gradient at velocity 
vd = cE c ff/B ~ vp. For graphene the drift velocity term 
in Eq. (1) can account for up to ~ 50% of the total and 
a measurement of EMP velocity offers an excellent probe 
of this important but hard to access quantity. 

The recent interest in bulk magnetoplasmons in 
graphene ji^-lHj has led to experiments on magnetic 
field splitting of the bulk plasmon frequency in the 
far-infrared absorption spectrum of spatially modulated 
graphene (2T 



25] 



The absorption resonance which was 
found to decrease in energy with increasing field is ex- 
pected to evolve into an edge-localised mode. However 
the propagation velocity, which is the central feature of 
EMPs, was not accessible in those experiments. 

In the present experiment, we create a wavepacket of 
EMPs with a voltage pulse and measure the time of prop- 
agation between two points on the edge such that left and 
right perimeter paths are different. Combining propa- 
gation time with path length measured by optical mi- 
croscopy affords a very direct measurement both of sign 
and amplitude of the group velocity which, from equation 
(I), is expected to be 
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Broadband microwave transmission was measured be- 
tween 2 points on the periphery of an oval shaped 
graphene flake of about f 5 x 5/zm exfoliated from natu- 
ral flake graphite onto the surface of a 290nm thick Si02 
layer commercially grown on an insulating Si wafer as 
shown in Fig. l(a,b). The flake is separated by ~ 3/im 
from the surrounding ground plane common to the two 50 
GHz bandwidth microwave coplanar waveguides whose 
2/im wide centre conductors overlap the edge by about 
2/im to serve as probe electrodes. One is capacitive and 
the other ohmic and they divide the perimeter in a ra- 
tio of about 2:1. The waveguides are tapered out to 
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FIG. 1: (a) Optical photograph of graphene sample with edge 
coupling electrodes terminating the coplanar waveguides, (b) 
Measurement configuration: EMP wavepacket excited by 7 or 
11 ps rise time 100 mV step function propagates along edge 
(white arrows), (c) Demodulated responses, offset for clarity, 
at B — ±18 T (filling factor v = 2) constructed by subtraction 
of waveforms at ±15V side gate potential for lips excitation. 
The arrival time difference arises from unequal left and right 
path lengths for oppositely directed magnetic fields. Overlay 
in red is best fit to weighted sum of first (blue dots) and 
second (solid blue) derivatives of a Fermi function (see text). 
Blue dot is at fitted arrival time: 21.4 ps for +18T and lips 
for -18T. (d) Direct measurement of 7 ps rise time excitation 
pulse broadened to 13 ps by the receiver amplifier. 



connect to a coplanar microwave printed circuit board 
placed horizontally and connectorised to vertical coax- 
ial lines leading to the room temperature excitation and 
signal analysis instruments. Total attenuation at 50GHz 
was -25dB. The sample was identified as graphene, and 
its quality monitored throughout fabrication, by scan- 
ning Raman spectroscopy (Supplementary Fig. SI). It 
was placed in cryogenic vacuum at 2.2K in a vertical 
magnetic field of up to ±19T where it was annealed and 
its density adjusted by heating to 425 K. Carrier den- 
sity was identified from Shubnikov-de Haas oscillations 
in microwave transmission at ~ 5 GHz, a technique anal- 
ogous to a 2-terminal DC measurement, but mobility 
was estimated from DC measurements on other, similarly 
prepared samples (26]. For the time of flight experiment 
pulse excitation and analysis are performed by a dual- 
channel time domain reflectometry module in transmis- 
sion mode. A fast 11 or 7 ps rise time step function pulse 
of ~ 100 mV (see Fig. 1(d)) is sent to one electrode 
to excite a wavepacket which propagates to the other 
electrode where it induces the signal sent to the receiver 
through a low noise 0.1-65 GHz preamplifier and digi- 
tally recorded by the picosecond sampling oscilloscope. 
The EMP wavepacket voltage is differentiated by the ca- 
pacitive coupling to give an approximately symmetric 
output pulse, but the raw waveform is a superposition 
of the response of the sample with the response of the 
transmission lines and amplifiers. To select the sample 



FIG. 2: (a) 3D colour plot of modulated signal amplitude in 
the propagation time - magnetic field plane, (b) Five signals 
chosen at filling factors ^ = 2,6, 10, offset for clarity. From top 
to bottom: 18T, 6T, 3.6T, 0T, -3.6T, -6T, -18T. These data 
correspond to lips excitation and 5ps sampling resolution. 



contribution, and to reduce noise from amplifier drift, a 
one hertz sidegate modulation voltage V sg — ±15V is ap- 
plied to the ohmic contact electrode via a broadband bias 
tee. The two signal waveforms are subtracted from one 
another to eliminate the response of all but the sample, 
it being the only element sensitive to V sg . The response 
was identical upon interchange of input and output lines 
with simultaneous reversal of magnetic field. 

The physics resides in the wave packet propagation 
times. Fig. 1(c) shows arrival signals for ±18T fields 
vs time measured from the leading edge of the exci- 
tation at the source module while Fig. 2 paints the 
overall picture of arrival time profiles for magnetic fields 
— 19 < B < +19 T corresponding to Landau level filling 
factors v — n s hc/eB > 1.9, where n s ~ 0.87 x 10 12 cm -2 
is the surface carrier density. Detailed analysis of the 
data proceeds by remarking that the side gate potential 
alters the filling factor at the edge, influencing both ve- 
locity and attenuation. Variation of attenuation mod- 
ulates amplitude and replicates the basic form of the 
wavepacket whereas variation of velocity modulates ar- 
rival time and differentiates the basic shape. The mean 
arrival times should be the same, motivating a fit to a 
weighted sum of symmetric and anti-symmetric parts de- 
rived from a single function. Guided by the exponential 
shoulders of the excitation step, we chose a Fermi func- 
tion for which a typical fit is illustrated in Fig. 1(c). Fig. 
3 plots arrival times deduced from the fits as a function 
of oppositely oriented magnetic fields. Field differentia- 
tion of the data calculated by subtracting arrival pulse 
waveforms at neighbouring fields corroborates the side 
gate modulated results (see Supplemental Material). 

The difference in arrival times for ± \B\ is compati- 
ble with the left and right edge path lengths of 14 and 
26 nm around the sample from emitter to receiver for 
fields above about 2.5 T. Converting arrival time into 
propagation time, however, requires knowing the time 



FIG. 3: Propagation times between emitter and receiver 
structures as function of perpendicular magnetic field for op- 
posing orientations. The determination of zero time is de- 
scribed in the text. Filling factors v = 10, 6 and 2 are at 3.6, 
6 and 18T respectively. Black (red) symbols correspond to 
the positive (negative) sign of the magnetic field. These data 
correspond to 7ps excitation and 2ps sampling resolution. 



of emission. We approach this problem in three ways. 
Firstly, the ratio of propagation lengths allows us to 
deduce time zero from the two arrival times supposing 
that the left and right speeds are the same at the same 
field intensity and that the wave propagates along the 
perimeter. Secondly, as B — > the propagation ve- 
locity v g oc a X y oc 1/B — > oo will be limited by the 
change in nature of the wave towards a zero field, non- 
chiral bulk plasmon, propagating with group velocity 

v g (B = 0) = yj!fv F ^ « ^vfv^v = 2), where q 
is of the order of the inverse width of the wave packet and 
kp = y/Wn~ s is the Fermi wavevector [28) . The arrival time 
is then the time required to travel at this much faster ve- 
locity along the direct internal path between emitter and 
receiver electrodes. Thirdly, if we neglect the formation 
of the bulk plasmon, the arrival times should extrapolate 
linearly to zero for B — >• 0. The three methods agree to 
±lps and fix the zero of the time scale on Fig. 3 from 
which we calculate the propagation velocities. 

Three significant features emerge from the propagation 
time data of Fig. 3. One is the positive, quantitative 
identification of chirality at high field, since the propaga- 
tion times are proportional to left and right path lengths 
for oppositely directed fields. Another is the unique pulse 
response at fixed field orientation indicating that only a 
single path is taken and no counter-propagating mode is 
detected. The third is that the propagation times and 
therefore velocities exhibit a plateau structure in mag- 
netic field with the same filling factor sequence as the 
Hall conductivity. Filling factors v = 10, 6 and 2 are es- 
timated to be at 3.6, 6 and 18T respectively. For the 
longer path the quantization is well developed at v = 2 
and is clearly discernable at v — 6. On the shorter path 
the quantization for v = 2 is observed, but for v = 6 it 
is somewhat blurred by the appearance of the dip fea- 
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ture between the two plateaus, the origin of which is not 
understood. 

The propagation velocities are determined by com- 
bining the propagation time results of Fig. 3 with the 
perimeter path lengths (as shown as function of the fill- 
ing factor in Figure 3 in the Supplemental Material). 
We extract the drift velocity vd{v) from the experimen- 
tal EMP propagation velocity at filling factors v = 2, 6 
and 10 using only the functional form in filling factor 
dependence of the two contributions - closely linear de- 
pendence for the Hall conductivity and the theoretically 
calculated ratios at mid gap chemical potential for the 



drift velo city [271 129(431|. The drift velocity so extracted 
is then independent of the fitting parameter qw. We find 
I'd = (0.7 ± 0.3) x 10 8 cm s _1 for v — 2, very close to 
the theoretical estimate vd — 0.7vf for v — 2 at a sharp 
edge. 

We estimate the product qw from the part of the prop- 
agation velocity linear in a xy by virtue of the logarithmic 
factor in equation ([2]) . To extract w from qw we estimate 
q by converting temporal width At given by the excita- 
tion rise time into spatial width As = v g Ar ~ 1/q. The 
received pulse widths are Ar p ~ 25ps of which about 20 
ps are due to the finite rise time of the excitation pulse 
combined with the receiver bandwidth. Replacing the 
lips with a 7ps excitation pulse made only a marginal 
difference of about 2ps to the received pulse width. Set- 
ting C — 1 in equation (2), strictly only appropriate for 
a step function edge profile, leads to w = 500 ± 200 nm, 
much greater than the expected £b ~ 6 nm. It could 
possibly be interpreted as lateral width of charge accu- 
mulation, but it could also result from not taking into 
account the x~ x l 2 density divergence on approach to the 
edge which, classically, accumulates charge (n ~ — v.Vn) 
of opposite sign to that of the edge itself lending a dipolar 
nature to the mode. 

EMP attenuation can be estimated from the ratio of 
response amplitudes for the two unequal paths if the at- 
tenuation coefficient is supposed the same for both. From 
the arrival signals for v = 2 shown in Fig. 1 the attenu- 
ation length is 70 ± 30 /jm corresponding to a relaxation 
time of 50 ± 20 ps, three orders of magnitude longer than 
the t ~ 0.05 ps Drude relaxation time applicable for bulk 
plasmons in similar samples pij which showed mobilities 
of- SOOOcm^V- 1 . 

The propagation times for the two paths are no longer 
experimentally distinguishable below about 2.5 T (v > 
14) , indicating that EMP excitations are no longer clearly 
chiral. From our experience on transport measurements 
in similar samples, these are fields for which full chiral- 
ity is no longer expected to be preserved as disorder su- 
presses the QHE gaps [26j|. Any echo from propagation 
around a full circle is at least 20 x smaller than the sin- 
gle path signal, indicating that at least one contact pad 
almost fully absorbs the wavepacket. For an ohmic con- 
tact, the low (50f2) impedance detector would effectively 



short-circuit the EMP wave impedance (p xy — 10 4 f2) 
and suppress re-emission. On the other hand, the second 
contact pad is capacitive because the main pulse has the 
shape of the derivative of the injected step voltage and 
there is no ohmic connection between the two edge elec- 
trodes. A model for chiral wave absorption by contact 
shunts is presented in Supplemental Material. 

In summary, graphene EMP differ from their inertial 
gapped 2DES counterparts not by the basic form of the 
dispersion relation but by the effect of the hard edge and 
gapless relativistic dynamics on the drift velocity. We 
identify them in graphene by observing chirality, quan- 
tization of velocity and low attenuation. We see no ev- 
idence of counter-propagating states under our experi- 
mental conditions, but these are far from the neutrality 
point for which such modes were proposed [32j]. The EMP 
propagation velocity exceeds the Fermi velocity, affirming 
the collective nature of the excitation. The carrier drift 
component vd is less, but comparable, to the Fermi veloc- 
ity, in keeping with gapless, massless particle behaviour 
at a sharp edge. The fit to a characteristic length of 500 
nm in the edge Coulomb interaction is one order of mag- 
nitude smaller than observed in soft edged GaAs[13] but 
it is nearly two orders of magnitude greater than the mag- 
netic length. Although this could reflect lateral charge 
distribution, it could also result from the oversimplified 
model which does not take into account the divergence 
of the edge density profile or electrostatic screening by 
the compressible strip. 

For use in plasmonics, bulk plasmons in graphene offer 
great advantages over metals in being electrically tun- 
able and in having long relaxation times (33l - [37| . Edge 
magneto-plasmons offer further features over bulk plas- 
mons: they are > 1000 x less damped, they are chiral 
with chirality determined by gate voltage sign and the 
propagation velocity is linearly rather fourth root depen- 
dent on gate voltage. Graphene EMP could extend con- 
siderably the frequency domain of graphene based plas- 
monics and open up a whole new domain of chiral plas- 
monics. 

In the course of completing the work reported here, we 
became aware of similar work on EMP propagation in 
mm size samples of graphene on SiC involving ns rather 
than ps propagation times |38j]. 
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Supplemental material: Carrier drift velocity and edge magnetoplasmons in 

graphene 



Sample Fabrication 



The sample was made by exfoliating graphene from natural Indian graphite procured from NGS Naturgraphit 
GmbH. It was deposited by the standard Scotch tape procedure on a 290nm thick layer of thermal SiC>2 commercially 
grown on an insulating Si wafer of room temperature resistivity ~ 8kr2-cm. The wafer was prepared in a Piranha 
solution (H2SO4 + H2O2) followed by an oxygen plasma etch. Spurious grains of graphite were removed with high 
pressure oxygen plasma etching, during which the graphene was protected with a PMMA/HSQ mask. The coplanar 
waveguide (CPWG) conductors of Ti(5nm)/Au(200nm) were then deposited by evaporation in high vacuum through 
an electron beam patterned PMMA mask. The central conductors extend onto the edge of the graphene by about 
2/im as seen in the photo of Fig. 1(a) of the main text. After liftoff the graphene sample was characterised with 
scanning micro- Raman spectroscopy (see Figure [4]) with a 532nm wavelength laser and found to be a monolayer with 
low impurity content. The wafer was cut into a 3 x 3mm 2 flip chip so as to leave the sample centrally placed and for 
the on chip CPWG to make contact with the CPWG of the microwave printed circuit board (PCB). 




FIG. 4: Raman sample characterisation, (a) Raman map of the sample, where green corresponds to graphene. (b) Raman 
spectrum in the middle of the graphene sample. 



Experimental configuration 

The sample was mounted perpendicular to the (vertical) field axis in a hermetic insert made for the bore of the 19T 
superconducting magnet. The magnet was cooled in a 4 He cryostat with lambda plate cooling of the coil region to 
~ 2.2K. The sample was annealed by heating its holder to 150°C for two hours while maintaining a cryogenic vacuum 
with the cold 4.2K walls of the insert. Transmission measurements of 52 1 as a function of field at a few GHz allowed 
us to estimate the electronic density to be about 0.87 x 10 12 cm~ 2 . 



Design of the microwave circuit 



We designed the microwave circuit to and from the sample to have 50GHz bandwidth with the help of CST 
Microwave Studio software. Two coaxial cryogenic cables connect two V-connector hermetic feed-throughs from the 
top of the insert cryostat to the horizontal sample holder PCB assembly via Mini-SMP right angle connectors modified 
to have less dielectric to ensure reduced reflection and with a bent central conductor to make solderless connections 
to the tapered coplanar waveguides (CPWG). The waveguides are etched onto a 17 /mi thick gold-on-copper face of 
a TMM10 substrate PCB with a dense array of via holes. The waveguides terminate at contact pads on adjacent 
edges of a 2 x 2mm 2 square hole over which is placed the 3 x 3mm 2 chip from an undoped Si wafer containing the 
graphene sample. We used undoped Si to limit dielectric losses. The "flip chip" , so called because it is mounted face 
down on the PCB to ensure contact between its own CPWG and those of the PCB, contains sample specific tapered 
waveguides which terminate at the graphene in either an ohmic or capacitive contact. The flip-chip is held in place 
by an optically flat metallic backing plate attached to a bent beam spring of Vespel plastic. All contacts to the PCB 
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are solderless to resist the high 450-500K temperatures required for annealing, care being taken to ensure sufficient 
elasticity of the contacts to compensate thermal contraction. The frequency-dependent loss in the PCB and flip-chip 
CPWG was measured by the transmission coefficient S21 between the two Mini-SMP connectors with a continuous 
CPWG on the flipchip. It has a linear baseline with small oscillations in frequency, and is about -7dB at 50GHz. The 
insertion loss of the total circuit including 4m of cryogenic coaxial cable is around 25dB at 50GHz. The complete 
circuit is operational up to 50GHz. At higher frequency vertical standing wave modes start to appear across the 
thickness of the PCB, leading to substantial loss. The details of the microwave setup will be published elsewhere. 



Magnetic field derivative of the transmitted signal 



60 
40 
20 


-20 



-15-10-5 5 10 15 

field (Tesla) 

FIG. 5: Field derivative of the transmitted signal. Rainbow colour plot of the derivative of the transmitted signal with 
respect to the magnetic field. Positive to negative values range across the spectrum from red to blue. The black and red 
superimposed data points are the transmitted pulse arrival times as shown in Figure 3 of the main text, where the signal is 
modulated with respect to gate voltage. 

The derivative with respect to field of the arrival signal is illustrated by the rainbow colour plot in Figure [5] which 
is obtained by subtracting consecutive 0.33T interval direct (unmodulated) transmission signals. The propagation 
times deduced from the analysis of the derivative with respect to gate voltage (Figure 3 in the main text) have 
been superimposed to facilitate comparison from which we conclude that the derivative in field gives similar pulse 
propagation times to those obtained from the derivative with respect to gate voltage. However, unlike the gate 
voltage modulated signal, the field derivative is more subject to experimental drift as the difference is made between 
acquisitions separated by long times and shows inferior signal to noise. 




Edge magnetoplasmon propagation velocity 
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FIG. 6: Propagation velocity as function of the filling factor. Propagation velocity as function of the filling factor for 
both signs of the magnetic field from the data presented on main text Figure 3 with the same colour coding. 
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Modeling the ac response of the graphene sample 

Below we derive the ac response of a graphene sample in the idealized case of perfectly transmitting edge channels 
and assuming that, at a given filling factor, all the chiral edge channels running in parallel propagate with the same 
velocity. The aim is to understand the effect of a perfect ohmic contact connected to an external circuit of negligible 
impedance on the chiral EMP propagation. Does it fully absorb the EMP chiral wave? How is the pulse response 
modified when the magnetic field is reversed? 

To model the ac response, we start from the gauge invariant, current conserving approach developed in [H which 
is particularly appropriate for chiral quantum Hall edge channels. We first derive the ac conductance for two simple 
situations: the case of two ohmic contacts and the case of one ohmic and one capacitive contact. The results of [l[ 
are valid in the low frequency limit where the long range interaction can be simply taken into account by introducing 
capacitances to ground. In the second part, we extend our analysis to include the high frequency regime relevant to 
edge magneto-plasmons by adapting the approach of Q for a ID wire with short range interactions to the chiral case. 
The assumption of chirality leads to a dispersion relation for the magneto-plasmon mode where the electron drift 
velocity simply adds to the plasmon velocity with short range interactions. Generalization to long range interactions 
would lead directly to Eq.(l) of the main text. Regarding the ac current response to a voltage pulse, our calculation 
shows that the only relevant modification due to interactions is to change the propagation velocity, while the absorptive 
condition of the contact remains unchanged. In particular the existence of at least one ohmic contact is shown to 
prevent the formation of pulse echoes and would seem responsible for our observing only a single pulse in the current 
response. 



/. low frequency regime 



contact 





2 


ih) 











contact 



FIG. 7: Two ohmic contacts. Conduction paths between two ohmic contacts connected by two chiral states of lengths h 
and h- 

We first consider the case of two ohmic contacts as shown in Figure [7J The lengths of the filling factor v chiral 
edge states connecting contacts (1) and (2) are l\ and 1%. In the limit where one can linearize the electron dispersion 
relation around the chemical potential of the edges, the ac transport Kirchhoff's laws relating the ac current I a at 
the input of contact a to the ac voltage Vp applied on contact /3 are 



h = v—[V 1 -V 2 e* 



V A {1 



(3) 



h = 



p 2 / ~ 
v— (-Vie 1 * 1 +V 2 -V 3 (1- e^ 1 ) 



(4) 



2 



(5) 
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2 

h = y\ (-V 2 (l - e^)+Vi(l - e*)) (6) 

where (f>x(2) — ^h(2)/ v D — l i- LUT i(2) is the temporal phase, uj the frequency, and vjj the carrier drift velocity. In 
the above set of equations the terms e 1 ^ 1 have to be replaced by 1 + iuiTi — l/2(u!Ti) 2 in the limit of small to. In 
this limit, the voltage V 3 (4j represents a fictitious voltage drop / 3 ( 4 )/(— icjC ?3(4) ) across the quantum capacitance 
^93(4) = e2rr i{2)/h- Upon including the geometrical capacitance C 3 ( 4 ), it is related to the actual voltage source by 

^3(4) = ^3(4) + ^3(4)/(~* w C3(4))- I n keeping with [l| equations (1-4) are invariant under a global shift of all potentials 
by a common quantity U and they show current conservation as they include the displacement currents 7 3 and I4. 

In the absence of backscattering between upper and lower edges, the dynamics of each edge are decoupled. If the 
two contacts are not equally spaced along the perimeter, as is intentionally the case for our graphene sample, i.e. 
li ^ l 2 , the ac response of the current is expected to change under magnetic field reversal. As an example, we take 
V\ = V while Vq^i = and look at the current response / = —h- This would correspond experimentally to the 
situation where the input transmission line is connected to contact (1) and the output transmission line to contact (2) 
and the finite characteristic impedance of the transmission line Zc = 50S1 is neglected with respect to h/ve 2 (several 
kil). To first order in uj 

Reversing the sign of the magnetic field, and therefore the direction of propagation, leads to a similar expression 
but with index (3) replaced by the index (4). Thus, unless the upper and lower lengths are equal, 1(B) ^ I(—B) 
Because of the displacement currents, the B parity imposed by the Onsager-Casimir relations for two-terminal 
samples no longer holds at finite frequency. 

In the limit of infinite geometrical capacitances V 3 (4) = U 3 (4) , the Coulomb interaction is fully screened and the 
following expression for I(V) is valid at any frequency 

/ = v e \e 1 ^ (8) 
h 

or in the time domain I(t) = v^-V(t — t±). If the magnetic field is reversed t% and <j)\ are replaced by t-i and 4>2- 




FIG. 8: One ohmic and one capacitive contact. Top and bottom figures differ by a gauge transformation which shifts 
the overall potential by U. 

We now consider the case where contact (2) is so resistive that it can be replaced by a small but finite capacitance 
Co connected to the output of the transmission line (Figure [S]). The non- interacting limit is treated by setting 

the geometrical capacitances C 3 and C4 to be infinite. The very short capacitive contact length (few fim) can be 

2 

represented by a quantum capacitance C q = v^-t in series with Co, small enough to be in the regime where lot <C 1 
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(r ~ few ps). Following 0, 0, Q the problem of the current 1 = —1% response to an ac voltage V can be replaced by 
an equivalent problem where the internal potential U in the sample is, according to gauge invariance, subtracted off 
all potentials, as in the lower part of Figure [51 Let l£ = vj^{V — U)e 1 ^' 1 be the current at the input of the capacitive 
contact, and Iq = I^e 1 ^ the current at the output of the contact. <j>i — ujt\ is the temporal phase accumulated 
by the current emitted by contact (1) up to the entrance of the capacitive contact and (f> = ujt the temporal phase 
accumulated within the capacitive contact. We have 

h = -I = / + (l - e*) = (-iC w) U. (9) 

For ujt < 1, we find the impedance 

T = + ^ (R* + -4rr) ■ (10) 

I —iCoUJ \ —iCqUJ J 

For 4>i — we recover the expression of the mesoscopic capacitor impedance of 0, H| , where R q = ■ For small ui 
and large Co we have I(t) ~ C q dV(t — Ti)/dt. We expect the dominant part of the pulse response to be the derivative 
of the step function pulse emitted by the generator used in the TDR mode. 



II. Inclusion of short range interactions 




FIG. 9: Two ohmic contacts. Conductance between two ohmic contacts connected with two chiral states of lengths h and 
h. 



We now introduce Coulomb interactions in the short range limit. We start with the simple case of two perfect 
ohmic contacts just as above in Section I. The problem being chiral, in the absence of backscattering we expect 
the upper and lower edge dynamics to be fully decoupled. We follow the approach used in Q for modeling the ac 
response of a (non-chiral) ID interacting wire. This approach treats the wire, or Tomonaga-Luttinger liquid, not 
in the framework of the bosonization technique but in the self-consistent random phase approximation (RPA). The 
Coulomb interaction is taken to be short-range. This is realized by assuming a capacitor placed very close to the 
electrons, see Figure [9] The interaction becomes local with ac potential u(x) = (e/c)p(x), where p(x) is the local 
charge density and c is the capacitance per unit length defined by the wire and the screening gate. The difference 
with respect to the non-chiral case enters in the polarization kernel H(x,x'). For the upper edge, the last right hand 
term of Eq. (4) in [2|, lqF £ F e nF\x-x ^ j g re placed by iqpVFe lqF{ - x ~ x ) for x' < x and by for x' > x as a result 
of chirality. Here we use the notation qp — lu/vd where the drift velocity replaces the Fermi velocity used for the 
non-chiral case. Vp = c q /e 2 is the one-dimensional density of states with c q = ve 2 /KVp) the quantum capacitance per 
unit length of the upper edge. Similar definitions hold for the lower edge. The solution of the self consistent equations 
for the propagation of the potential or the charge density along the upper edge is sought in the form of e tqx (while 
both components e ±lqx are necessarily present in the non-chiral case). We find q — qp-^r. As a result, the phase 
velocity of the, now collective, excitation describing the charge propagation becomes 

a ve 2 ( 1 1 \ 

W0 = — = —r- \- - ) =v D + vemp (11) 

q h \c q cj 
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We recover the expression of Eq.(l) in the main text where the edge magneto-plasmon velocity in the screened 
regime is now Vemp — Qxy/c- 

Coming back to the problem of the current at contact (2) in response to the ac voltage applied at contact (1), the 
only modification to the expression obtained in Section I is to replace the drift velocity by the velocity v^. Specifically, 
in the notation of Figure [71 we take V2 — and V4 — 0, such that the two edges are decoupled. We find 



We see that this is the same solution as the one given by equations (1-3) with <j>\ replaced by ql\. Adding the 
contribution of a finite V2 and V4, a set of equations identical to Eq.(l-4) is recovered with V$u} replaced by VW) 
(indeed, the geometrical capacitances are already taken into account by the short range interaction which renormalizes 
the propagation velocity). The main consequence for the pulse response is that the role of the ohmic contact is the 
same as in the non-interacting case. Despite the strong change in the dynamics of the current propagation, full 
absorption of the pulse is expected if at least one ohmic contact is present. This prevents the pulses from making 
several turns along the sample edge and no echo is expected, as is presently the case in our experiment. 
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